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A  Note  on  the  Brewster  Angle  in  Lossy  Dielectric  Media 


In  electromagnetism,  the  Brewster  angle,  dB,  is  the  angle  of  incidence  at  which  there  is  no 
reflection  of  the  parallel  (i.e.  vertical)  polarized  wave  from  a  planar  interface  between  materials 
with  differing  permittivity.  That  is,  the  reflection  coefficient  goes  to  zero  at  this  angle  and  all  of 
the  incident  power  is  transmitted  into  the  dielectric  medium.  The  reflection  coefficient  is  a 
function  of  the  constitutive  parameters  permittivity,  e,  and  permeability,  /v,  of  the  two  media 
forming  the  interface  as  well  as  the  angle  of  incidence.  For  media  with  /urt  =  /un ,  the  Brewster 
angle  is  commonly  expressed  as 


where  the  r  subscript  denotes  the  relative  constitutive  parameter,  the  /  subscript  denotes  the 
parameter  of  the  medium  in  which  the  incident  wave  is  propagating,  and  the  t  subscript  denotes 
parameter  of  the  medium  into  which  the  wave  is  transmitted.  We  assume  the  incident  medium 
is  free  space  and  the  transmitted  medium  is  soil  for  the  purposes  of  ground  penetrating  radar 
(GPR).  For  perpendicularly  polarized  waves,  no  real  incident  angle  exists  that  will  reduce  the 
reflection  coefficient  to  zero.  [1] 

For  clarity,  we  now  define  some  basic  definitions  associated  with  plane  waves  incident  at  a 
planar  interface: 

1.  Plane  of  incidence  -  the  plane  containing  the  propagation  vector  of  the  incident  wave  /c, 
and  the  unit  normal  to  the  interface. 

2.  Parallel  polarization  -  the  electric  field  of  the  incident  wave  lies  in  the  plane  of  incidence. 
The  orthogonal  magnetic  field  is  in  the  direction  transverse  to  plane  of  incidence,  so  we 
call  this  a  transverse  magnetic  (TM)  wave.  This  is  referred  to  as  vertical  polarization  in 
the  GPR  community. 

3.  Perpendicular  polarization  -  the  electric  field  of  the  incident  wave  lies  normal  to  the 
plane  of  incidence.  Since  the  electric  field  is  in  the  direction  transverse  to  the  plane  of 
incidence,  we  call  this  a  transverse  electric  (TE)  wave.  This  is  referred  to  as  horizontal 
polarization  in  the  GPR  community. 

In  realistic  environments  involving  GPR,  one  must  consider  the  lossy  component  of  the  soil’s 
permittivity.  Soils  become  lossy  as  moisture  is  added,  especially  in  mineralized  soils  with  high 
salinity  [2],  [3].  The  soil  moisture  strongly  affects  the  real  part  of  the  permittivity.  While  the 
salinity  of  soils  does  not  have  a  dramatic  influence  on  the  real  part  of  the  permittivity,  the 
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imaginary  part  is  strongly  affected  by  both  the  salinity  and  moisture  of  the  soil,  especially  at  low 
frequencies. 

Expressing  the  complex  permittivity  for  the  incident  wave  medium  and  the  transmitted  wave 
medium  as 


(2) 

(3) 


we  can  attempt  to  find  the  Brewster  angle  using  (1)  as  follows.  We  change  the  complex 
permittivity  into  phasor  form  to  make  the  mathematics  associated  with  the  square  root  slightly 
easier.  That  is,  we  can  write 


0B  =  tan  1  — — — —  =  tan  1 

\e'--  is" 

\  ri  J  n 


<p(.M, ) 


/  ~  r  -i 

=  tan^  -^expy-(4-4)  . 


where  we  use  cp  to  represent  the  phase  angle  as  opposed  to  the  azimuthal  angle  of  incidence. 
Now,  changing  back  from  phasor  form  in  order  to  take  the  inverse  tangent,  the  argument 
becomes 


We  then  attempt  to  solve  for  the  Brewster  angle  using  the  definition  of  the  inverse  tangent  of  a 
complex  number  as 


0B  =  tan  '(a  +  jb ) 


zJwM 

2  \\  +  j{a  +  jb)\ 


(not  a  realizable  angle!)  (6) 


where  the  constants  a  and  b  are  defined  in  (5)  and  In  is  the  natural  logarithm. 


It  can  be  seen  from  (6)  that  the  Brewster  angle  for  a  complex  valued  permittivity  is  a  complex 
value.  Since  complex  valued  physical  angles  do  not  exist,  we  can  conclude  that  the  reflection 
coefficient  cannot  truly  go  to  zero  in  a  soil  with  complex  valued  permittivity.  In  other  words,  in  a 
lossy  soil  there  is  no  angle  for  which  all  of  the  incident  TM  wave  power  will  be  transmitted  into 
the  soil. 
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Though  the  reflection  coefficient  does  not  exactly  go  to  zero  in  a  lossy  medium,  there  still  exists 
a  minimum  value  at  some  angle,  which  is  often  referred  to  as  the  pseudo-Brewster  angle.  We 
define  the  reflection  coefficient  for  a  TM  wave  incident  upon  a  soil  from  free  space  as 


ru  n 2  sin  9, .  -  Jn 2  -  cos 2  9 , 

T™  =  — - '  ^  ,  (7) 

n 2  sin  9 1  +-sJn2  -  cos2  9l 

where  0,  is  the  incident  angle  and  n  is  the  complex  index  of  refraction  of  the  soil 

n  =  V(4  -je”rt'htr,)  ■  (8) 

If  we  assume  a  nonmagnetic  medium,  the  relative  permeability  in  (8)  will  be  unity,  i.e.  prf  =  1 . 
We  can  observe  how  the  imaginary  component  of  the  permittivity  affects  the  wave  propagation 
at  a  hypothetical  interface  by  plotting  the  reflection  coefficient  for  several  values  of  imaginary 
permittivity  while  keeping  the  real  part  constant.  Figure  1  shows  the  magnitude  of  the  reflection 
coefficient  plotted  versus  angle  of  incidence  for  four  different  imaginary  permittivity  components. 
The  real  component  of  the  permittivity  was  kept  constant  at  unity.  For  the  case  of  no  loss  and  a 
very  low  (unity)  real  permittivity,  we  observe  no  reflection.  For  the  lossy  case,  we  can  see  that 
the  imaginary  component  of  the  permittivity  affects  the  angle  at  which  we  find  the  dip  in  the 
reflection  coefficient.  We  can  also  observe  that  this  dip  becomes  sharper  for  higher  losses.  We 
can  draw  two  important  conclusions  from  these  observations: 


1 .  One  cannot  simply  use  the  real  part  of  the  permittivity  to  determine  the  angle  for  minimum 
reflection  (which  we  may  call  the  “pseudo-Brewster  angle”)  given  a  lossy  soil. 

2.  In  real  environments  the  soil  surface  (dielectric  interface)  is  not  perfectly  flat.  If  there  is  a 
large  amount  of  loss  in  the  soil,  small  deviations  in  angle  around  the  pseudo-Brewster  angle 
will  significantly  increase  the  reflection  coefficient. 


As  an  example  demonstrating  the  second  point  stated  above,  the  change  in  the  reflection 
coefficient  due  to  angle  deviations  around  the  pseudo-Brewster  angle  can  be  observed  in  Figure 
1  for  the  highest  loss  case.  The  minimum  value  of  the  reflection  coefficient  for  the  case  where 

an  =  1  -  ylOO  occurs  when  r  =0.17at  0,  =  84.25°.  The  reflection  coefficient  doubles  if  the 

angle  of  incidence  is  reduced  by  -8.75°  or  if  it  is  increased  by  only  3.5°.  This  change  in 
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reflection  due  to  small  deviations  in  angle  would  preclude  the  exploitation  of  the  Brewster  angle 
for  GPR  applications  in  lossy  soils  with  realistic  undulating  surfaces. 

However,  if  we  examine  the  case  of  a  medium  with  low  loss,  say?rt  =1  -j,  the  reflection 

coefficient  doubles  from  its  minimum  if  the  angle  is  reduced  by  -16.85°  or  if  it  is  increased  by 
10.05°.  In  the  low  loss  case,  there  is  slightly  more  leeway  for  angle  deviations  while  keeping  the 
reflection  coefficient  relatively  low. 


Figure  1 .  Reflection  coefficient  for  media  of  varying  levels  of  loss. 


When  moisture  is  added  to  soil,  the  real  part  of  the  permittivity  is  affected  as  well  as  the 
imaginary  component.  Using  data  collected  from  the  U.S.  Army  RDECOM  CERDEC  Night 
Vision  and  Electronic  Sensors  Directorate’s  mine  lanes  facility  and  reported  in  [4],  the  reflection 
coefficient  can  be  determined  for  a  soil  with  varying  levels  of  moisture.  The  soil  under  test  in  the 
present  analysis  is  called  bluestone  gravel.  The  complex  permittivity  was  derived  using  S- 
parameter  measurements  of  a  soil  sample  in  a  coaxial  holder  from  1  MHz  -  2  GHz.  Varying 
levels  of  moisture  were  added  to  the  sample  and  the  volumetric  moisture  percentage  was 
recorded.  Table  1  shows  the  permittivity  values  reported  from  the  measurements  at  200  MHz. 
Figure  2  shows  the  reflection  coefficient  for  the  bluestone  gravel  versus  incident  angle  as 
calculated  using  equation  (7)  for  a  frequency  of  200  MHz. 
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Table  1.  Permittivity  values  of  “bluestone  gravel”  soil  from  the  NVESD  mine  lanes  facility  at  200  MHz. 


Measured  Permittivity  for  Bluestone  Gravel  at  200  MHz 

Volumetric  Moisture  (%)  Permittivity  (srt  =  s'rt  -  js"n ) 

0.8 

3.0-  j0.2 

2.5 

3.6  -  j0.7 

8.3 

5.1  -  j2.7 

15.8 

00 

1 

'oo 

o 

30.3 

19.0  —  j  1 7.0 

Figure  2.  Reflection  coefficient  for  soil  from  the  NVESD  mine  lanes 
facility  for  varying  levels  of  moisture  at  200  MHz. 
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Table  2  shows  the  permittivity  values  reported  from  measurements  at  2  GHz.  Figure  3  shows 
the  reflection  coefficient  for  the  bluestone  gravel  versus  incident  angle  as  calculated  using 
equation  (7)  for  a  frequency  of  2  GHz.  While  the  reflection  coefficient  reaches  a  minimum  very 
close  to  zero  in  this  case,  it  should  be  noted  that  this  minimum  changes  drastically  with  changes 
in  volumetric  moisture. 


Table  2.  Permittivity  values  of  “bluestone  gravel”  soil  from  the  NVESD  mine  lanes  facility  at  2  GHz. 


Measured  Permittivity  for  Bluestone  Gravel  at  2  GHz 

Volumetric  Moisture  (%)  Permittivity  ( srt  =  s'n  -  js"rt ) 

0.8 

2.8-j0.0 

2.5 

3.1  -jO.1 

8.3 

4.4  -  j0.5 

15.8 

7.5  —  j  1 .4 

30.3 

17.5 -j3.3 
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Figure  3.  Reflection  coefficient  for  soil  from  the  NVESD  mine  lanes 
facility  for  varying  levels  of  moisture  at  2  GHz. 

The  Brewster  angle  expression  given  in  (6)  merely  shows  that  there  is  no  physical  angle  at 
which  the  reflectivity  goes  to  zero  for  a  lossy  medium  and  should  not  be  used  for  actual 
calculations.  One  may  be  tempted  to  take  the  real  part  or  the  absolute  value  of  the  complex 
angle,  but  neither  of  these  corresponds  to  the  minimum  value  of  the  reflection  coefficient.  It  can 
be  shown  that  the  reflection  coefficient  reaches  a  minimum  when  the  real  part  of  r  =  0 .  As 
stated  in  [5]  without  derivation,  the  pseudo-Brewster  angle,  which  is  the  angle  for  which  the 
absolute  magnitude  of  the  complex  reflection  coefficient  reaches  its  minimum,  0B  ,  is  given  as 


for  s'ri  =  land  s"ri  =  Oas  in  the  case  where  the  incident  field  is  in  free  space.  The  angle  given  in 

(9)  is  the  angle  from  grazing.  Equation  (9)  can  be  verified  by  examining  the  reflection 
coefficients  plotted  in  Figure  2.  For  example,  the  calculated  minimum  reflection  coefficient  at 

200  MHz  for  30.3%  volumetric  moisture  (s  =19-/17)  is  Irl"  =78.74°,  which  is  in 

accordance  with  the  plotted  value.  The  calculated  minimum  reflection  coefficient  at  2  GHz  for 

30.3%  volumetric  moisture  (art  =17.5-y3.3)  is  |F| ^  =76.66° .  This  example  not  only  shows 

the  validity  of  equation  (9),  it  also  shows  that  dispersion  plays  a  role  in  determining  the  pseudo- 
Brewster  angle.  Since  the  complex  permittivity  is  a  function  of  frequency,  the  reflection 
coefficient  must  be  a  function  of  frequency  too.  Most  GPR  systems  are  wideband,  which  means 
that  the  pseudo-Brewster  angle  will  not  be  constant  over  the  entire  radiated  band. 

It  has  been  shown  in  this  note  that  the  Brewster  angle  does  not  exist  for  dielectric  interfaces 
involving  complex  permittivity.  A  pseudo-Brewster  angle  does  exist,  which  is  the  angle  at  which 
the  magnitude  of  the  reflection  coefficient  reaches  a  minimum.  However,  the  pseudo-Brewster 
angle  cannot  be  calculated  using  the  standard  Brewster  angle  formula.  Furthermore,  the 
reflection  coefficient  can  rapidly  rise  for  small  angle  deviations  away  from  the  pseudo-Brewster 
angle  when  the  soil  is  very  lossy.  Finally,  the  dispersive  nature  of  lossy  soils  means  that  the 
pseudo-Brewster  angle  will  not  be  constant  over  a  broad  band. 


A  Note  on  the  Brewster  Angle  in  Lossy  Dielectric  Media;  Ian  McMichael,  US  Army  RDECOM  CERDEC  NVESD 

7 


References 

[1]  C.  Balanis,  Advanced  Engineering  Electromagnetics,  New  York:  John  Wiley  &  Sons,  1989. 

[2]  Y.  Shao,  Q.  Hu,  H  Gao,  Y.  Lu,  Q.  Dong,  C.  Han,  “Effect  of  Dielectric  Properties  of  Moist 
Salinized  Soils  on  Backscattering  Coefficients  Extracted  From  RADARSAT  Image,”  IEEE  Trans. 
Geosci.  Remote  Sensing ,  vol.  41,  no.  8,  pp.  1879-1888,  Aug.  2003. 

[3]  J.  Hipp,  “Soil  Electromagnetic  Parameters  as  Functions  of  Frequency,  Soil  Density,  and  Soil 
Moisture,”  Proc.  of  the  IEEE,  vol.  62,  pp.  98-103,  Jan.  1974. 

[4]  J.  Curtis,  D.  Leavell,  C.  Weiss,  R.  North,  E.  Smith,  J.  Cortes,  R.  Castellane,  M.  Fields, 
“Characterization  of  Soils  from  the  Night  Vision  and  Electronic  Sensors  Directorate  Mine  Lane 
Facility,  Fort  Belvoir,  VA,”  U.S.  Army  Corps  of  Engineers,  Engineer  Research  and  Development 
Center,  July  2003. 

[5]  E.  Shotland,  “Three  Angles  Significant  in  Radio  Propagation,”  IEEE  Trans.  Antennas 
Propag.,  vol.  20,  no.  6,  pp.  798-801,  Nov.  1972. 


A  Note  on  the  Brewster  Angle  in  Lossy  Dielectric  Media;  Ian  McMichael,  US  Army  RDECOM  CERDEC  NVESD 

S 


Appendix 

Derivation  for  the  pseudo-Brewster  angle  for  media  with  complex  valued  permittivity: 


The  reflection  coefficient  for  a  parallel  polarized  wave  incident  from  free  space  onto  a 
nonmagnetic  medium  is  given  as 


where  n  =  Js'rt  -  js”rt  for  nonmagnetic  media.  With  complex  valued  permittivity  and 
nonmagnetic  media,  we  can  rewrite  (A1)  as 


Y™  — 


(K  -  X  )sin  -  Vfe  -  je'r, ) -  cos2 

(K,  - JK )sin 0,  +  J(s'rl  -j s"rt)-cos2di 


Each  square  root  can  be  separated  into  a  complex  number  using  the  following  identity: 


where  the  sign  should  be  chosen  to  be  the  same  as  the  sign  of  y.  After  separating  the  square 
roots,  we  can  combine  real  and  imaginary  terms  and  then  multiply  the  top  and  bottom  by  the 
conjugate  of  the  denominator  to  get  the  following 


The  pseudo-Brewster  angle  for  lossy  media  0B ,  which  we  are  defining  as  the  angle  for  which 
the  magnitude  of  the  complex  reflection  coefficient  reaches  its  minimum,  occurs  when  the  real 
part  of  the  complex  reflection  coefficient  equals  zero.  Therefore,  to  solve  for  the  Brewster 
angle,  we  simply  set  the  real  part  of  (A4)  equal  to  zero.  The  denominator  of  the  real  part  of  (A4) 
vanishes  when  solving  the  algebraic  equation  so  that  we  only  need  to  solve  the  following 
equation  for  0, 


t'r 2  +  Kt1  )sin2  dt  -  V(4  -  COS 2  e)2  +s"rt  =  0  .  (A7) 

Using  the  double  angle  formula  and  combining  like  cosine  terms,  (A7)  takes  the  form 
cos2  267 

a  b  c 


k2+<2J 


-1 


+  cos  2  0, 


4s'.,  -  2 (e’rt~  +  s"rt 2 )  -  2  +  (s'rt2  +  s"rt 2 )  -  4 (s', 2  +  s"2 )+  4 s'rt  - 1 


=  0 


(A8) 


We  then  use  the  quadratic  equation 


cos  267  = 


—  b  ± 

2  a 


-  4  ac 


(A9) 


where  the  constants  are  defined  as  shown  in  equation  (A8).  We  choose  the  root  with  the 
negative  sign  in  the  numerator  to  assure  that  the  quotient  is  between  -1  and  1  and  corresponds 
to  real  incident  angles.  The  final  step  is  to  take  the  inverse  cosine  and  solve  for  0„  which 

satisfies  the  condition  foxdB  ,  as  follows 


On  =  —  COS  1 


-b-^Jb2  -4 ac 


2  a 


(A10) 


The  expression  given  in  (A10)  can  be  condensed  by  assigning 

E  =  s'2  +  sf  , 

which  simplifies  the  pseudo-Brewster  angle  to 


(All) 


9n  =  —  cos  1 


(l  +  £2  -2s') 

i-V 

(2s: -E2 

->) 

r- 

(e2  -\\e2  -4E  +  4s'r-\) 

1 

(e2 

->) 

(A12) 
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The  solution  (All)  can  be  shown  to  be  equivalent  to  that  of  (9),  which  was  given  in  the 
literature. 
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